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1. Introduction 

Non-perturbative Hamiltonian light-front quantum field theory presents opportunities and chal- 
lenges that bridge particle physics and nuclear physics. Major goals include predicting both the 
masses and transitions rates of the hadrons and their structures as seen in high-momentum trans- 
fer experiments. Current focii of intense experimental and theoretical research that could benefit 
from insights derived within this Hamiltonian approach include the spin structure of the proton, the 
neutron electromagnetic form factor, the generalized parton distributions of the baryons, etc. 

Hamiltonian light-front field theory in a discretized momentum basis ^ and in transverse 
lattice approaches [Q, |3|] have shown significant promise. We present here a basis-function approach 
that exploits recent advances in solving the non-relativistic strongly interacting nuclear many-body 
problem [Q, H]. We note that both fields of physics face common issues within the Hamiltonian 
approach - i.e. how to (1) define the Hamiltonian; (2) renormalize to a finite space; (3) solve for 
non-perturbative observables while preserving all symmetries; and, (4) take the continuum limit. 

We begin with a brief overview of recent advances in solving light nuclei with realistic nucleon- 
nucleon (NN) and three-nucleon (NNN) interactions using ab initio no-core methods. Then, we 
introduce our basis function approach to light-front QCD within the light-front gauge. Renormal- 
ization/regularization issues are also addressed. We present illustrative features of our approach 
with the example of cavity-mode QED and sketch its extension to cavity-mode QCD. 

2. No Core Shell Model (NCSM) and No Core Full Configuration (NCFC) methods 

To solve for the properties of nuclei, self-bound strongly interacting systems, with realistic 
Hamiltonians, one faces immense analytical and computational challenges. Recently, ab initio 
approaches have been developed that preserve all the underlying symmetries and they converge to 
the exact result. The basis function approach that we adopt here ||] is one of several methods 
shown to be successful. The primary advantages are its flexibility for choosing the Hamiltonian, 
the method of renormalization/regularization and the basis space. These advantages impel us to 
adopt the basis function approach in light-front quantum field theory. 

While non-relativistic applications in finite nuclei restrict the basis to a fixed number of fermions, 
we introduce here the extension to a flexible number of fermions, antifermions and bosons. 

Refs. [Q] and ^ provide examples of the recent advances in the ab initio NCSM and NCFC, 
respectively. The former adopts a basis-space renormalization method and applies it to realistic 
NN and NNN interactions (derived from chiral effective field theory) to solve nuclei with Atomic 
Numbers 10-13 Experimental binding energies, spectra, electromagnetic moments and transi- 
tion rates are well-reproduced. The latter adopts a realistic NN interaction that is sufficiently soft 
that binding energies from a sequence of finite matrix solutions may be extrapolated to the infinite 
matrix limit. Again, one obtains good agreement with experiment. 

It is important to note the analytical and technical advances made to solve these problems. 
First, non-perturbative renormalization has been developed to accompany these basis-space meth- 
ods and their success is impressive. Several schemes have emerged and current research focuses 
on detailed understanding of the scheme-dependence of convergence rates (different observables 
converge at different rates) [f7|]. Second, large scale calculations are performed on leadership-class 
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parallel computers to solve for the low-lying eigenstates and eigenvectors as well as to carry out 
evaluation of a suite of experimental observables. Low-lying solutions for matrices of dimension 
2-billion on 30,628 processors with a 3-hour run is the current record. However, one expects these 
limits to be exceeded very soon as the techniques are evolving rapidly [||] and the computers are 
growing dramatically. Matrices with dimensions in the tens of billions with strong interaction 
Hamiltonians will soon be solvable. 

In a NCSM or NCFC application, one adopts a 3-D harmonic oscillator for all the particles in 
the nucleus (with harmonic oscillator energy HQ.), treats the neutrons and protons independently, 
and generates a many-fermion basis space that includes the lowest oscillator configurations as well 
as all those generated by allowing up to N max oscillator quanta of excitations. The single particle 
states specify the orbital angular momentum projection and the basis is referred to as the m-scheme 
basis. For the NCSM one also selects a renormalization scheme linked to the basis truncation 
while in the NCFC the renormalization is either absent or of a type that retains the infinite matrix 
problem. In the NCFC case [g], one extrapolates to the continuum limit as we now illustrate. 
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Figure 1: Calculated ground state energy of 12 C for N max = 2—8 (discrete points) at selected values of HQ.. 
For each h£l, the results are fit to an exponential plus a constant, the asymptote, which is constrained to be 
the same for each curve|Q|. We display the experimental ground state energy and the common asymptote. 

We show in fig. |l| results for the ground state of C as a function of N max obtained with 
a realistic NN interaction, JISP16 [ph. The smooth curves portray fits that achieve the desired 
independence of N max and HQ. so as to yield the extrapolated ground state energy. Our assessed 
uncertainty in the extrapolant is about 2 MeV so there is rather good agreement with experiment. 
The largest cases correspond to N max = 8, where the matrix reaches a dimension near 600 million. 
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N max = 10 produces a matrix near 8 billion which we plan to solve in the coming year. 



3. Cavity mode light-front QED and QCD 



It has long been known that light-front Hamiltonian quantum field theory has similarities with 
non-relativistic quantum many-body theory. We further exploit this connection, in what we will 
term a "Basis Light Front Quantized (BLFQ)" approach, by adopting a light-front basis space 
consisting of the 2-D harmonic oscillator for the transverse modes (radial coordinate p and polar 
angle 0) and a discretized momentum space basis for the longitudinal modes with anti-periodic 
boundary conditions (APBC). The 2-D oscillator states are characterized by their principal quantum 
number n and orbital quantum number m. Adoption of this basis is also consistent with recent 



developments in AdS/CFT correspondence with QCD []1 OfJ - 

In order to bring the full 3-D basis into perspective, we select a transverse and the lowest 
longitudinal mode (with index j = 0) and display slices of the 3-D basis function at selected lon- 
gitudinal coordinates, x~ in fig. ||. Our purpose in presenting figs. ^ is to suggest the richness and 
economy of texture available for solutions in a basis function approach. Note that the choice of 
basis functions is arbitrary except for the standard conditions of orthonormality and completeness. 



APBC: k + =-(j + -) 
L 2 




- L < x < +L 
j = 0, n = 1, m = 




Figure 2: Transverse sections of a 3-D basis function involving a 2-D harmonic oscillator and a longitudinal 
mode with antiperiodic boundary conditions (APBC). 

To further illustrate the BLFQ approach, we consider a non-interacting QED system confined 
to a transverse harmonic trap or cavity. For simplicity, we take the leptons as massless. The 
basis functions are matched to the trap so we implement a transverse 2-D harmonic oscillator 
basis with length scale fixed by the trap and an infinite system in the longitudinal direction with 
APBC. The symmetries (and their chosen values for the case illustrated in fig. |3|) imposed on the 
many-parton basis states are total charge (Z=3), total magnetic projection (M = 0) in the transverse 
space, total spin projection (S = 1/2), maximum total 2-D oscillator quanta (N max = 8 — 19) and 
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total longitudinal momentum in dimensionless units j, + 1 = K = 6). The range of the number 
fermion-antifermion pairs and bosons is limited by the cutoffs in the modes (N max and K): e.g. each 
partem carries at least one unit of longitudinal momentum so the basis consists of up to 6 partons. 

Fig. f| demonstrates the saturation of low-lying unperturbed modes with increasing N max . 
Already, this state density could serve as input to a model for the statistical mechanics of the system 
treated in the microcanonical ensemble. Of course, interactions must be added to make the model 
realistic at low temperatures where correlations are important. After turning on the interactions, the 
challenge will be to evaluate observables and demonstrate convergence with respect to the cutoffs 
{Nmax an d K). Independence of the basis scale, Ml, must also be obtained. These are the standard 
challenges of taking the continuum limit. 

The Hamiltonian H for this system (KH gives the invariant mass-squared) is defined by the 
sum of the occupied modes with the scale set by the combined constant A 2 = TMqWL: 

H = 2M P- = — - — V 1 1 
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Figure 3: State density from BLFQ for 3 identical massless fermions confined in a trap for a selection of 
Nmax values at fixed K = 6. The dimensions of the resulting matrices are presented in the legend. The states 
are binned in groups of 5 units of energy (quanta) where each parton carries quanta equal to its 2-D oscillator 
energy divided by its light-front momentum fraction (x, = (j; + 1)/K). The dashed straight line ttaces the 
exponential increase in state density, familiar in many-body theory, when the basis is sufficiently large. 

We can extend the approach to QCD by implementing the SU(3) color degree of freedom for 
each parton - 3 colors for each fermion and 8 for each boson. We consider two versions of imple- 



menting the global color singlet constraint. In the first case, we follow Ref . [11] by constraining 
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all color components to have zero color projection and adding a Lagrange multiplier term to the 
Hamiltonian to select global color singlet eigenstates. This results in the upper curves in fig. f|. 
In the second case, we restrict the basis space to global color singlets and this results in the lower 
curves in fig. [|. The second method produces a factor of 30-40 lower multiplicity at the upper ends 
of these curves at the cost of increased computation time for matrix elements. Either implemen- 
tation dramatically increases the state density over the case of QED, but the use of a global color 
singlet constraint is clearly more effective in minimizing the explosion in basis space states. 
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Figure 4: Number of color space states that apply to each space-spin configuration of selected multi-parton 
states for two methods of enumerating the color basis states. The upper curves are counts of all color 
configurations with zero color projection. The lower curves are counts of global color singlets. 



Now we briefly address the interacting theory. We have identified the appearance of the ex- 
pected log divergence in the matrix element of the fermion-boson vertex as it falls only as the 
inverse square root of the principal quantum number in selected sequences. We are investigating 



sector-dependent coupling constant renormalization to manage this divergence [12]. 



4. Conclusion 



Following successful methods of ab initio nuclear many-body theory, we have introduced a 
basis light-front quantization (BLFQ) approach to Hamiltonian quantum field theory and illustrated 
some of its key features with a cavity mode treatment of QED. We indicated our method of treating 
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color for extending our approach to QCD. The computational requirements of this approach are 
substantial, and we foresee extensive use of leadership-class computers to obtain practical results. 
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